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Abstract 

The c-axis quasiparticle conductivity shows different behavior depending on 
the nature of the interlayer coupling. For coherent coupling with a constant 
hopping amplitude t±, the conductivity at zero frequency and zero tempera- 
ture <r(0, 0) depends on the direction of the magnetic field, but it does not for 
angle-dependent hopping t((f>) which removes the contribution of the nodal 
quasiparticles. For incoherent coupling, the conductivity is also independent 
of field direction and changes only when paramagnetic effects are included. 
The conductivity sum rule can be used to determine the admixture of coherent 
to incoherent coupling. The value of <r(0, 0) can be dominated by t± while at 
the same time t{4>) dominates the temperature dependence of the superfluid 
density. 

PACS numbers: 74.20.-z, 74.25. Gz 
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I. INTRODUCTION 



The nature of the interlayer coupling between two adjacent Cu0 2 planes in the cuprates is 
an important issue that remains unresolved. Suggestions include effects of strong intralayer 
scattering,!! non- Fermi liquid ground states,! the general phenomenon of confinement,! inter- 
and in-plane charge fluctuations,! indirect c-axis coupling through the particle-particle 
channel,! as well as resonant tunneling on localized states in the blocking layer! and two 
band models.0 Coherent coupling originates from an overlap of the electronic wave functions 
between planes, and in-plane momentum is conserved in interlayer hopping. By contrast 
for impurity-mediated incoherent coupling, the in-plane momentum is not constraint. It 
has been shown that the c-axis conductivity sum rule depends on the nature of interlayer 
coupling.! Coherent coupling obeys the conventional sum rule! regardless of the angular 
dependence of the interlayer hopping amplitude. On the other hand, incoherent coupling 
violates the sum rule even if the in-plane dynamics can be described by a Fermi liquid. 
However, in order to explain the violation observed in some experiments]^ it was also nec- 
essary to include the non-Fermi liquid nature of the in-plane dynamics.0 For YBa2Cu307_,5 
at optimum doping, a conventional c-axis sum rule is observed^ which is consistent with 
coherent c-axis coupling and an in-plane Fermi liquid. For the underdoped case, a pseudogap 
is observed and the sum rule is closer to l/2.0lll This value can most easily be understood 
as a pseudogap effect with incoherent c-axis coupling.0 

Another interesting example of the interlayer coupling is the c-axis conductivity due to 
low lying quasiparticles at zero temperature and zero frequency. A method for observing the 
quasiparticle current is to measure the hysteretic c-axis tunneling current- volt age charac- 
teristics of layered cuprates.0 For measurements of the optical conductivity which includes 
the pair tunneling, see Ref.0 For coherent coupling with a constant hopping amplitude t±, 
the c-axis quasiparticle conductivity a{w — > 0,T — > 0) shows a residual value indepen- 
dent of the in-plane scattering rate 7 C A (gap) to leading orderll! with the correction 
— (7/A ) 2 . The underlying physics in this case is the same as for the in-plane quasiparticle 
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conduct ivity;lljil3 namely, the impurity- induced density of states at uj = is canceled by 
the decrease of quasiparticle lifetime. A universal value has been observed in YE^CusOg.g 
for the in-plane thermal conduct ivity.S In some sense the thermal conductivity is an ideal 
probe of universal behavior since, as opposed to the electrical conductivity, it is not renor- 
malized by vertex and by Fermi-liquid corrections,^ and so there is less ambiguity in its 
identification. The universal conductivity limit does not appear for coherent coupling with 
an angle- dependent hopping amplitude of the form t(4>) = t0cos 2 (20) which is believed 
to be appropriate for the copper oxides, where <fi = tan -1 (k y /k x ) is an angle in the mo- 
mentum space. In this case the c-axis quasiparticle conductivity is reduced by a factor of 
(3/8) (7/A0) 2 as compared with the value of a(w — > 0, T — > 0) for a constant t± assuming the 
same magnitude of hopping amplitude, which of course is not the case. This arises because 
the angular dependence, cos 2 (20), eliminates the contribution of the quasiparticles on the 
nodal lines from interlayer transport. For incoherent coupling, the residual conductivity is 
proportional to (7 m[A /7]/A ) in leading order. Consequently, the universal value of the 
c-axis quasiparticle conductivity is a characteristic only of coherent coupling with a constant 
hopping amplitude. Geometrical consideration of the Cu and O atomic arrangement s^B 
from plane to plane leads one to expect that the dominant overlapping of orbitals would lead 
to a t(4>) form with possibly a subdominant constant piece t±. In this case, t± would still 
dominate the value of cx(0, 0) for very small 7, but the other terms would become important 
as (7/A0) 2 corrections become larger in which case the term t(<j>) or the incoherent part can 
become important. 

In this paper we investigate effects of interlayer couplings on the c-axis quasiparticle 
transport in the absence, as well as in the presence of an in-plane magnetic field. It is 
our aim in this paper to look at all three contributions of interlayer coupling in detail, and 
our discussion of the magnetic field effects will not be restricted to the constant t± as it is 
111 Ref.0 We also investigate the role played by the conductivity sum rule in quasiparticle 
transport. We find that it can be used to fix the relative amount of coherent to incoherent 
coupling. We also find that ct(uj, 0)/a(0, 0) shows a different uj 2 coefficient for uj < 7 and 



cr(0, T)/er(0, 0) a different T 2 coefficient (T < 7) depending on the nature of the coupling. 
In the presence of an in-plane magnetic field, cx q (0, 0) for coherent coupling with a constant 
hopping amplitude t± depends on the direction of the field0 as does the coefficient of the 
T 2 term in <r q (0,T), but such a dependence is negligible for an angle-dependent coherent 
t(4>) because it removes the contribution from the nodal quasiparticles which otherwise 
manifest the effect of field direction. For incoherent coupling, <r q (0, 0) is independent of field 
direction, and in fact changes only when a paramagnetic interaction is included. Here we 
mention that there exist, in the literature, studies which consider effects of oxigen doping on 
the c-axis transport and Coulomb charging effect! associated with a pseudogap behavior.il 
Such complication, however, is beyond the scope of the present paper. 

This paper is organized as follows: In Sec. II, we derive general formulas associated with 
the c-axis quasiparticle conductivity including paramagnetic effects. We describe, in Sec. Ill, 
that effects of coherent coupling with a constant as well as with an angle- dependent hopping 
amplitude on the c-axis quasiparticle transport with and without an in-plane magnetic field. 
In Sec. IV, impurity-mediated incoherent c-axis coupling is considered. We also illustrate the 
role of the conductivity sum rule in quasiparticle transport in Sec. V and draw conclusions 
in Sec. VI. 



II. FORMALISM 

The Hamiltonian H for a cuprate superconductor with interlayer coupling can be written 
as H = H + H c , where H describes a d-wave superconductor and 

# c = £ [t k - P c+ (k)c 2CT ( P ) + h.c] . (1) 

cr,k,p 

The interlayer hopping amplitude tk- P depends on nature of c-axis couplings: i) tk- P = 
£±$k-p for coherent coupling with a constant amplitude, ii) For an angluar dependent am- 
plitude tk-p = £<^k-p cos 2 (2(f)). In the lattice, it can be seen from geometrical consideration^ 
that £^<5k-p [cos(k x a) — cos(A^a)] 2 with an in-plane lattice constant a. ii) For incoherent cou- 
pling tk-p = Vk-p and an impurity average is to be taken into account. 



Applying perturbation theory for H c , we obtain the c-axis quasiparticle conductivity 
a(u,T) in terms of the electronic spectral weight function A(k,e). 

a{u) = -- E *k- P / rfe W + - /(e)] ^(k, e) A(p, 6 + W ) , (2) 

k,p 

where C is a constant which depends on the nature of interlayer coupling, /(e) is the usual 
Fermi thermal factor, and the spectral function 

A(k A - 7(1 + Wgjj 7(1 - Wgk) , , 



with = Y^k + ^k- We assume Ak = A [cos(k x a) — cos(k y a)} (or A cos(2</>) in the 
continuum limit) and a cylindrical Fermi surface with £k = k 2 /(2m) — e F measured from 
the Fermi energy e F . As T ->• and w -»• 0, <x(0,0) = C E k , P *k- P ^( k > °)^(p> 0). It is 
easy to calculate cx(0,0)(= cr ) for the couplings we mentioned earlier. To obtain the c-axis 
quasiparticle conductivity at finite frequency to at zero temperature, we take T — > limit 
and assume u; < 7 



k,p 



<U(0-i (£)#/(«0 + yj^ 3 /( e 



2\ 

' A(k,e)A(p,e + u) 



6=0 



^CY: tl_ p A(k, 0)A(p, u)-C£ 4-p (%) d.A(k, e)A(p, e + u) 
k, P k, P vz/ 

+ C E 4-p f ^ ) d 2 £ A(k, e)A(p, e + u;) + • • • (4) 



where an expansion in uj has been made. The finite T corrections at u = will be considered 
later along with an in-plane magnetic field. 

In the presence of an in-plane magnetic field, we assume that the field penetrates freely 
into the sample so that the field is uniform between the planes.0 This assumption should be 
valid for Bi- and Tl-based cuprates which are nearly two dimensional. For YBCO, however, 
there is a possibility that in-plane vortices form. This means that in such a case an average 
over vortices in a unit cell is required.il The interlayer coupling Hamiltonian is modified by 
the presence of a uniform in-plane magnetic field H as follows: 



H c = £ [t k - P C+ (k + q)C 2CT (p) + h.c] , (5) 



where q = (erf/2) (z x H) with an interlayer distance rf. 

Following a procedure similar to the one we applied in the zero field case, we obtain the 
c-axis quasiparticle conductivity cr q (u;,T) in the presence of an in-plane magnetic field, 

^) = -~ £ £-p / & [/(c + w) - /(e)] A(k + q, e)A(p, e + u) . (6) 

k,p 

In our consideration, the energy scale associated with the magnetic field are always much less 
than the gap A . Due to an in-plane field, the quasiparticles gain an additional momentum 
q on transferring to the next plane. In other words, the Bogoliubov-de Gennes (BdG) wave 
functions^ Uk(r) and t>k(r) become ■ui c e*( k+q - ) ' r and t>ke l( - k+q )' r , respectively.il The spectral 
weight function becomes 

1 J (c - £ k+q ) + 7 2 (e + £ k+q ) 2 + 7 2 1 ^ 



where -Ek+q — y £k+ q + ^k with £k+q — (k • q)/m. Note that the angular dependence 
of Oq comes from the factor (k ■ q)/m ~ E q cos(<f) — 9 q ), where E q = Vpq and 9 q is the 
direction of q, which can be interpreted as the direction of the field (9) because of the 
symmetry in the problem. The gap seen by quasiparticles with a momentum k + q is not 
Ak but Ak+ q ; however, Ak+ q = Ak + 5 q , where <5 q = qv g cos(0 + 9 q ) with vq — 2A /&f, 
and the correction <5 q makes a negligible contribution to the quasiparticle energy spectrum. 
When the paramagnetic interaction is included, -Ek+q — > Ek+ q ± HbH in the denominator 
of A(k + q, e), where \xb is the Bohr magneton. As to — > and T —>■ 0, 

cx q (cu - 0, T - 0) = C £ 4_ p A(k + q, 0) A(p, 0) . (8) 

k,p 

Based on symmetry consideration,^ one can easily deduce that the dependence of 
<T q (u; — > 0,T — > 0) = cr(#) on field angle has a period 7r/2 because of the rf-wave gap. 
The sign of the gap does not matter. Consequently, 

a{9) = Cn cos(4n9) , (9) 

n=0 

where the C n 's depend on the ratio of the magnetic field energy to 7 i.e., on Eg/j and hbH/j. 
Since both E q and fisH are linear to H, we can parameterize the ratio (E q / hbH) ~ d/a. We 

6 



choose {E q l HbH) = 6 for the high-T c Bi- and Tl-based cuprates.il Another relation can be 
deduced without detailed calculation. If we expand a (9) in terms of E q /^, then its angular 
dependence will appear for the first time from the term (Eg/j) 4 cos 4 (</> — 9) because the 
fourth harmonic cos (46*) comes from the cos 4 (#) term. This means that 7 < E q is required 
for any clear angular dependence to show up. For the eighth harmonic cos(8#) to enter, the 
field has to be even higher. 

At a finite temperature T < 7, we apply the Sommerfeld expansion' to er q (u; — > 0, T) = 
a(9, T) and obtain 

. T) = -C £ ti_ p J de (^fj A(k + q, e)A(p, 



a(0, T) = -C £ f£_ B / de ( ) A(k + q, e)A(p, e) 

,p 



= a(0) 



2 / m\ 2 



6 I 7 



a{9) 



(10) 



where 



grEic,pt£- p ^ + q,^(p,O e< „ 

" l J E k ,ptL P ^(k + q,0)A(p,0) ' 1 J 

with e' = e/7. 

We use the nodal approximation" to describe the low temperature physics associated 
with the quasiparticle transport because quasiparticles near the nodal lines dominantly 
contribute to the resistive transport. Linearizing £k and Ak near the nodal points on the 
Fermi surface (FS), we obtain £k = f^(k • kj_ — hp) and Ak = u^k • kn, where vp(kp) is the 
Fermi velocity (momentum) and v G = \Z2aA is the gap velocity.' The unit vector (kn) 
is perpendicular (parallel) to the FS. These unit vectors will change depending on the nodal 
line. For example, on the nodal line of <p = 7r/4, kj_ = (k^+kj,) / y/2 and ky = (— k^+kj,) / a/2. 
Since we include quasiparticles only near the nodal points on the FS in our considerations, 
we will use the following procedure in the actual calculation. 

E-E/^-E/mm*. (12) 

k node v ' node 



where p\ = vpk_L, P2 = vaku and J = \{2t[) 2 vfVg\ 1 ■ The coordinate transformations 
we made are rotation, translation and dilation. In the coordinate of {pi,P2), the energy 



dispersion of the quasiparticle E^ becomes Jp\ + p\ < po, where p ~ O(A ). 

III. COHERENT COUPLING 

In this section, we consider the effects of coherent coupling on c-axis quasiparticle trans- 
port. Since tk- P = t±5k- p for a constant hopping amplitude, as uj — > and T — > 0, 
°o = CEk^l^dc, o) 2 . For uj < 7, we use the nodal approximation for Eq. (fj). Then, up to 
order (cj/7) 2 , 

— = 1 + U-)'< < 13 > 

a 18 \ jj 

where o"o is a constant independent of 7. 

In the presence of a uniform in-plane field, the quasiparticle conductivity becomes 

a{6) = C^t 2 L A(k + q,0)A(k,0) . (14) 

k 



Applying the substitution Eq. flT2| ) to the above equation, we obtain without the paramag- 
netic part 

pdp 
p 2 + 1 

(15) 



a(0) = t\CJ £ J dip J 
node 

1 



p 2 + 2(E 9 /7) cos(0 - e)pcos(if) + (^g/7) 2 cos 2 (0 - 9) + 1 ' 

It can be shown within the nodal approximation that = <p n ± tan" 1 [^/^g^f], where <p n 
is a direction on the nodal line; however, to a good approximation, we can replace <fi by <p n . 
For the weak field case E q < 7, we expand a (6) in terms of E q /^y and obtain 



^ 12 V 7 / 80 



1 - - cos(4#) 



. 7 



1 



(16) 



where <To is the c-axis quasiparticle conductivity in the absence of the in-plane magnetic 
field. As we see, the angular dependence of a (6) is small for a weak field. a (6) is maximum 
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(minmum) when the in-plane field is along a nodal (anti-nodal) line. It can be physically 
understood as follows: When H is along a nodal line, for example, 9 = 7r/4, the angle of q 
is 37r/4 because q oc (z x H). Then, while the shifted momenta k + q of the quasiparticles 
with cf) = 7r/4 and <p — 5tt/4 deviate from the nodal regions when q is compatible with k, 
quasiparticles with = 37r/4 and <fi = 7ir/4 remain in the nodal areas. This means the 
remaining quasiparticles govern the c-axis transport when 6 = 7r/4. On the other hand, for 
6 = 0, all quasiparticles move away from the nodal regions due to the momentum shift q. 
Therefore, a{6) < <7o and cr(0) < c(7r/4) because of a mismatch between k and k + q in the 
interlayer transport. 

In Fig. 1, we plot resistivity as a function of direction of the in-plane magnetic field, 
p{6) = l/a(6) for various values of E q with (solid line) and without (dashed line) the 
paramagnetic interaction included. We reproduced results of Ref.0 for E q /^f = 1, 2, and 
4. As shown in Fig. 1, p(8) is increased as the magnetic field increases, and it is decreased 
for a give E q when the paramagnetic interaction is considered because the interaction is 
pair breaking. If we consider only the paramagnetic interaction, then the quasiparticle 
conductivity a has no 9 dependence and increases as follows: 



The eighth harmonic clearly appears for E q /^ = 12. Actually, it begins to appear when 
Eq/j ^ 10. Interestingly, the paramagnetic interaction unambiguously reduces the ampli- 
tude of p(9) for a higher field (E q 3> 7). We also plot a{9) of Eq. (pi]) for a finite T < 7 in 
Fig. 2. For reference in the absence of the magnetic field, a(8) =4/3 (dashed line). 

In the case of an angle- dependent hopping amplitude tk- P = £(0)<5k- P , the c-axis quasi- 
particle conductivity in zero field is a Q = C2 k t(^) 2 A(k, 0) 2 . For uj < 7, using the nodal 
approximation we obtain 



where o"o oc (7/A0) 2 in this case. Effects of the in-plane field on the c-axis quasiparticle 
conductivity can be seen in the same way as before, namely, 




(17) 




(18) 
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a(6) = CY, t(<P) 2 A(k + q, 0)A(k, 0) . (19) 

k 

Since t(<p) = ^(j^/Ao) 2 within the nodal approximation, a{6) becomes without the param- 
agnetic interaction 

pdp ( p \ 4 



node 

(20) 



p 2 + 2(^/7) cos(0 - 9)pcos(<p) + {E q /-i) 2 cos 2 (0 - 6) + 1 ' 

In this case, we find that cr(8)/ao — 1 without clear dependence on the direction of the 
in-plane field, and a (9) ~ O [(7/A0) 2 ln(A /7)] in the leading order and the next order is 

[(E q / Aq) 2 }; therefore, a(9) is negligible. The c-axis quasiparticle conductivity is insen- 
sitive to the in-plane magnetic field because the angle- dependent hopping amplitude t±((p) 
eliminates contributions of nodal quasiparticles, which otherwise manifest the effects of the 
field. 

IV. INCOHERENT COUPLING 

For impurity- mediated incoherent coupling tk-p = Vk-p, we need a model for the impu- 
rity scattering potential and need to carry out an average over impurity configurations. We 
use a simple modeli! for the scattering potential |Vk_ p | 2 = | Vo | 2 + |Vi| 2 cos(20 fe ) cos(20 p ). 
One may expand |Vk- p | 2 with respect to scattering symmetry so that it is decomposed 

|^k-p| 2 = I Vol 2 + £ \V lv \ 2 cos[(4/ + 2)<f> h ] cos[(4/' + 2)</> p ] , (21) 
i,v 

where I and I' are integers. However, for the c-axis quasiparticle transport, only the 

1 Vo| 2 term gives a non-zero value to the conductivity. As uj — > and T — > 0, cr = 
C\V \ 2 Ek, P ^( k , 0)^4(p, 0). For u < 7, it can be shown that 

°M 1 1 1 (22) 
<7„ 31n(A / 7 )' 1 ' 

where ctq oc [7ln(A /7)/A ] 2 for impurity- mediated incoherent coupling. 
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In the presence of an in-plane magnetic field, the c-axis quasiparticle conductivity be- 
comes 



= ciHI 2 E ^( k + q. °M(p> °) • ( 23 ) 

k,p 

As we see, the angular dependence of a{9) is determined by Sk^4( k + q>0). It can be 
shown numerically and analytically that the field effect appears only as the paramagnetic 
contribution. We present some analytic results, for example, in the weak field limit 

E ^( k + q 0) - ^( k > 0)1 * h (^-) , (24) 

where /i(x) = x 2 — x 4 /6 + a; 6 /15. Therefore, the c-axis quasiparticle conductivity for a weak 
field becomes 

a{9) ^ h^ B Hh) 

v. ~ + ln(A / 7 ) • {2b) 

Consequently, for incoherent coupling the in-plane magnetic field has no effect on the c-axis 
quasiparticle transport if the paramagnetic interaction is not considered. Physically, the 
momentum of the transport quasiparticle is not constraint in incoherent coupling so that 
the change in momentum q does not matter. Since the paramagnetic interaction is pair 
breaking, the magnitude of a (9) is increased but shows no angular dependence when this 
interaction is included. At a finite temperature (T < 7), we use Eq. (pi]) and for a weak 
field, a{9) becomes 

a (Q) ~ I + t^ 2 / 3 ) (26) 

{ ) + [ln(A / 7 ) + h{fi B H/i)) [1 - ^sHhf] ■ 1 ) 



V. C-AXIS CONDUCTIVITY SUM RULE 

From the c-axis conductivity sum rule, the superfluid density p s can be written in terms 
of the missing spectral weight (N n — N s ) and the theral averages of kinetic energy (H c ) s ^ 
of a superconducting (normal) state as follows: 
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p s = (N n - N s ) - Ane 2 d (H c ) s - (H c ) n 



(27) 



where u c is the cutoff frequency for the interband transitions that H c does not describe. 
Using the above equation as well as the Kramers- Kronig relation between the conductivity 
and the penetration depth, we obtain the normalized missing spectral weight (N n — N s ) / p s ; 
(N n - N s ) _ 1 lE 1 ,E k ,p|tk-p| 2 [G'(k,^)G'(p,^)-G' (k,^)G'o(p,^)] 



2 + 2 



(28) 



Ps * * £ w £ kiP |t k _ p | 2 F(k,^)F+(p,ccO 

where G(k,u>) and F(k,u) are superconducting Green functions and G (k, uS) is in the 
normal state.i0 

Since, in principle, all possible interlayer couplings might be present in a given sample, 
we need to consider coherent and incoherent coupling at the same time. The actual c-axis 
transport of quasiparticles is determined from competing effects between couplings. Based 
on the conductivity sum rule, we can estimate the ratio {o~ co /cr in ), where a co (t7j n ) is the 
c-axis quasiparticle conductivity for coherent (incoherent) coupling. Since the sum rule does 
not distinguish between t± and t(<f>), we consider two limiting cases: i) tk- P = £_L^k-p + ^k- P 
and ii) tk- P = £(0)<5k- P + ^k- P - As u> — > and T — > 0, for the case with t± 



a = o r 



Or 



a a, 



1 + 



= cr, 
4A a, 



1 + ^n 

Cm 



^ln(Ao/7)) : 



(29) 



ne 2 dt\N(0) \A 

where we have used a co = 2e 2 dt 2 ± N(0)/(7rA ) and a in = (8/7T 2 )a cn [7 ln(A /7)/A ] 2 with 
cr cn = 47rnj<i(eVoA^(0)) 2 . One may think that we need to know the ratio of |Vo| 2 to t\ in 
order to compare a co with <Ji n . However, the conductivity sum rule helps us to estimate the 
ratio ((Jco/ a in) without ad hoc information. 

Let us define 77 = p s , co / Ps,in = (iMDAV^iri)' where p s and A are the corresponding 
superfluid density and penetration depth, respectively. For coherent coupling 

1 1. 



1 

A 2 



l67Te 2 dt 2 ± N(0) 



-K(iA Q / 7 ) 



(30) 



where K is the complete elliptic integral of the first kind. For incoherent coupling, we assume 
\Vo 1 2 = |K| 2 for simplicity and to illustrate possible effects. Then, 
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^ = 32a m y £ [{k' 2 /k)K{k) - E(«)/« 



(31) 



where k = A /y A 2 , + (a; + 7sgnaj) 2 , k' = yl — k 2 , and E is the complete elliptic integral 
of the second kind.0 Since 7 << A , rj ~ 87re 2 c/t 2 L A^(0)/[12A a cn ] and we obtain 



a = a,, 



1 



(8/3) / 7 WA . ; 2 ' 



(32) 



Since the total superfTuid density p s = p SjCO + p Sj i n , it can be shown that 

N n -N s l + 2r] 1.08 



(33) 



p s 2(1 + 77) I + 77 

Note that r/ in Eq. (|32|) can be determined from the violation of the conductivity sum rule 
Eq. (|33|). rj —>■ 00 corresponds to pure coherent c-axis coupling in which case the sum rule 
Eq. (|33|) is conventional and equal to 1. On the other hand, r] — > means pure incoherent 
coupling and the sum rule is larger than 1 and equal to 1.58 in our model with |Vo/Vi| = 1. 
A finite t] corresponds to an admixture of coherent and incoherent c-axis transport. For 
example, rj ~ 5 applies when (N n — N s )/p s = 1.1, that is, there is a 10% violation of the 
sum rule upwards. The sum rule itself determines the admixture of p SjCO to p s ,im but cannot 
differentiate between t± and t(<p) for the coherent part. Limits on the relative size of these 
two overlap integrals can only be set from consideration of the chemistry of the CuC>2 planes 
and their overlap or, alternatively, from experimental information such as the behavior of 
cr(0, 0) when the impurity content is increased. 
For the case with t(<p) 



4- ^ 167re 2 ^ 2 ,AT(0) 



7 



=E 



A 



(34) 



16 S7T^2 + A 2 \^2 + Al 

where E is the complete elliptic integral of the second kind, and rj as above with 3 <-> 8 and 
t(j, <r->t±. Following a similar way as before, we obtain 

, (8/3) 



n 



MA0/7) 



(35) 



where a' co = <j co with (3/4) <-> 2 and t$ t±. 
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Note that t± and case both contribute in the same order to the zero temperature 
c-axis superfluid density. For the pure case, assuming t± and t$ are of the same magnitude, 
the ratio of the two contributions is 8/3. This is in sharp contrast to what we found in the 
case of <r(0, 0) where t± contributes of order t]_ while the contribution from t$ is of order 
/^(7/Ao) 2 and so vanishes in the clean limit. Recently, Gaifullin et a/.0 have measured 
the temperature dependence of the c-axis superfluid density in Bi2212 and found that at 
low temperature it is well fit by a form [1 — A(T/T c ) a ] with a of order 4 — 6 close to the 
values reported for Hgl20lS This finding favours a pure t((p) model which is known to 
give a T 5 law (Ref.@) and leaves little room for a subdominant t± contribution because 
in this case the c-axis penetration depth will mirror perfectly!! its in-plane temperature 
dependence which goes like T. The data certainly preclude a linear T contribution to the 
superfluid density of more than a few percent implying a ratio (t±/t ( j > ) of order < 1CT 1 at 
most. In this instance, the ratio (7/A0) can easily be comparable to (t±/t$) and can even be 
the larger of the two. As an illustration in the experimental work establishing the in-plane 
universal limit for the thermal conduct ivityS the values of (7/A0) are of order < 10 _1 . 
Latyshev et a/.0 quoted 7 ~ 3meV for their single crystals, which is of similar order. This 
implies that in the experiments on Bi2212, the contribution to er(0, 0) of each of t± and 
tj, terms are likely to be close in magnitude. Nevertheless, it is concluded in Ref.0 that 
the universal limit is observed so that t± presumably still dominates. Additional work with 
various values of 7, somewhat larger as well as smaller than used so far, should allow one to 
establish the size of the important ratio (t±/t^). This information, with the c-axis optical 
sum rule, would allow an estimate of each of the three possible contributions to the c-axis 
transport discussed in this paper. These parameters are important because they control 
all c-axis transport. We point out here that we have considered separately the case t± and 
t(4>). If amplitudes are added and squared, a cross term of the form 2tj_t^ cos 2 (20) will 
enter but this will not change qualitatively the conclusion made above. The contribution of 
such a term to a is (4e 2 GL/V(0)/7rA )t_i_t 9 i(7/Ao) 2 ln(A /7) and to the superfluid density is 
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167r 2 eW(0)t ± ^ 



1/2 - (2/7r)( 7 /^7 2 + A§)E( 7 /V7 2 + A§) 



VI. CONCLUSIONS 



We have discussed how the nature of the interlayer coupling influences the c-axis quasi- 
particle conductivity cr q (u;,T) in the absence and in the presence of an in-plane magnetic 
field. In zero field, cx(0, 0) is independent of the in-plane scattering rate 7, to leading or- 
der, only for coherent coupling with a constant hopping amplitude. Also, cr(u, 0) shows a 
different behavior for different c-axis coupling. Similarly, the field effects depend crucially 
on the nature of the interalyer coupling. For coherent coupling with a constant hopping 
amplitude, an angular dependence of the conductivity a(9) appears in high field. The re- 
sistivity p(9) increases because the mismatch between k and k + q decreases cr(9), and p(9) 
is maximum (minimum) when the field is along the anti-nodal (nodal) line. This confirms 
previous work.0 When the paramagnetic interaction is included, p{9) is decreased because 
this interaction is pair breaking and the anisotropy also decreases. For high field (E q ^> 7), 
the eighth harmonic appears and can be seen in the 9 dependence of p(9). For coherent 
coupling with an angle-dependent hopping amplitude, a (9) is insensitive to field direction 
because this hopping amplitute eliminates the contribution of the nodal quasiparticles which 
manifest the angular dependence. Also, the paramagnetic interaction does not significantly 
change u{9). For incoherent coupling, only the paramagnetic interaction has an effect on the 
c-axis quasiparticle transport and a(9) has no dependence on field direction. The coefficient 
of the first temperature correction which goes like T 2 is anisotropic only in the constant t± 
case and shows the same minima (maxima) as does the leading contribution. 

The c-axis conductivity sum rule helps estimate separately the contributions from coher- 
ent and incoherent coupling to the quasiparticle transport without microscopic information 
on hopping amplitudes. However, it cannot, by itself, differentiate between contributions 
from constant or angular dependent coherent hopping amplitude. To get separate infor- 
mation on these two contributions, impurity studies could be used. Consideration of the 



15 



temperature dependence of the c-axis superfluid density measured in Josephson plasma 
resonance experiments already gives evidence that is much larger than t±, but do not 
unambiguously rule out a small t± contribution which could still be large enough to dom- 
inate the universal limit for cr(0, 0, ) in relatively pure samples (t±/t ( j ) > 7/A0). At the 
same time, t& would dominate the temperature dependence of the superfluid density. We 



mention that the sum rule of Eq. ( p3|) was derived under the assumption that the normal 
state is Fermi-liquid like. However, if it has a non-Fermi-liquid nature such as a pseudogap, 
the sum rule has to be modified to account for this; therefore, it is necessary to explore the 
competing effects of interlayer coupling in a more fundamental theory which goes beyond a 
Fermi liquid normal state. 
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FIGURES 

FIG. 1. For various values of E q /j(= 1,4,6,12), resistivity p(9)/po = (Tq/u{9) as a function 
of a direction of the in-plane field 9 is plotted with (solid curve) and without (dashed curve) 
the paramagnetic interaction. As E q /^ is increased, p(9)/po is increased, and for a high field 
(Eg/ 7 = 12), the paramagnetic interaction unambiguously reduces the amplitude of p(9) and its 
anisotropy 



FIG. 2. a{9) in Eq. (10) shows finite temperature effects on the c-axis qusiparticle conductivity 
for Eg/'y = 1, 2, and 3. The paramagnetic interaction is included. In the absence of the field, 
a{8) = 4/3 (dahsed line). 
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Fig. 1 (Kim et al.) 




Fig. 2 (Kim et al) 




